Kitaev honeycomb and other exotic spin models with polar molecules 
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We show that ultracold polar molecules pinned in an optical lattice can be used to access a 
variety of exotic spin models, including the Kitaev honeycomb model. Treating each molecule as a 
rigid rotor, we use DC electric and microwave fields to define superpositions of rotational levels as 
effective spin degrees of freedom, while dipole-dipole interactions give rise to interactions between 
the spins. In particular, we show that, with sufficient microwave control, the interaction between 
two spins can be written as a sum of five independently controllable Hamiltonian terms proportional 
to the five rank-2 spherical harmonics Y2, q (6, </>), where (6, (j>) are the spherical coordinates of the 
vector connecting the two molecules. To demonstrate the potential of this approach beyond the 
simplest examples studied in [S. R. Manmana et al, arXiv:1210.5518v2], we focus on the realization 
of the Kitaev honeycomb model, which can support exotic non-Abelian anyonic excitations. We 
also discuss the possibility of generating spin Hamiltonians with arbitrary spin S, including those 
exhibiting SU(N=2S+1) symmetry 



INTRODUCTION 

Recent experimental progress in the control of ultra- 
cold polar molecules [1-6] has stimulated the study of 
many-body systems featuring strong dipole-dipole in- 
teractions [7-10]. A particularly promising research 
area [11 24] of implementing lattice Hamiltonians from 
dipolc-dipole-interacting molecules has emerged follow- 
ing the pioneering works of Rcfs. [25, 26]. This re- 
search area largely owes its promise to the great degree of 
controllability that external DC, microwave, and optical 
fields can provide over rotational levels and dipole-dipole 
interactions [11-15, 18, 26-32]. In Ref. [33], we demon- 
strated that this controllability can yield a variety of spin 
Hamiltonians, in which spin states are encoded in super- 
positions of rotational states of the molecules. While 
dipole-dipole interactions are often thought of as being 
proportional to 1 — 3 cos 2 8, where R = (R, 6, 4>) are the 
spherical coordinates of the vector connecting the two 
dipoles, we showed that not only the overall amplitude 
and sign of the spin Hamiltonian but also its individual 
terms can depend on R. We harnessed this dependence 
to propose the realization of spin models known to harbor 
topological phases of matter [33]. 

In contrast to Ref. [33] , which focused on the simplest 
examples, in the present paper we discuss a range of mod- 
els that are more interesting but harder to implement ex- 
perimentally. Specifically, as in Ref. [33], we emphasize 
that the interaction Hamiltonian between two molecules 
can, in principle, be a sum of five independently con- 
trollable terms proportional to the five rank-2 spherical 
harmonics Y2 iq (0,cf)). 

We then study a specific class of Hamiltonians with an 
arbitrary spin S. This class of Hamiltonians can be used 
to realize [33] the general SU(2) symmetric 5=1 model, 
i.e. the bilinear-biquadratic spin model, which has a rich 



phase diagram even in one dimension [11, 34, 35]. It can 
also give rise to SU(A^=2S , +1) spin models, which have 
recently attracted a great deal of theoretical and exper- 
imental attention due to their potential to access exotic 
strongly-correlated phases in alkaline-earth atoms [36- 
39]. The strong dipole-dipole interactions of the polar- 
molecule implementation of SU(A^) models can be larger 
than the SU(A f )-symmetric superexchange interactions 
in cold-atom implementations [36]. 

Then we turn to the primary example of this paper: 
the implementation of the Kitaev honeycomb model. In 
contrast to the proposal of Ref. [28] , our implementation 
of the Kitaev honeycomb model relies on direct - rather 
than perturbative - dipole-dipole interaction and hence 
gives rise to stronger interactions, which are easier to 
access experimentally. Other proposals for implementing 
the Kitaev honeycomb model use solid-state Mott insula- 
tors [40]. trapped ions [41]. superconducting quantum cir- 
cuits [42], weak superexchange interaction between cold 
atoms in optical lattices [43] , and an experimentally chal- 
lenging system of coupled cavity arrays [44] . We end the 
manuscript with a discussion of how we use microwave 
fields to create the dressed states that enable our imple- 
mentations, followed by an outlook. 

SETUP 

Consider a deep optical lattice in the XY plane shown 
in Fig. 1 and loaded with one polar molecule per site. 
Each molecule is described by a rigid rotor with rota- 
tional angular momentum operator N, rotational con- 
stant B, and dipole moment operator d. In the presence 
of an electric field E along z, each molecule is described 
by the Hamiltonian 

H = BN 2 - Ed z . (1) 
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FIG. 1. (color online). The lattice of molecules is in the 
XY plane. The direction of the DC electric field is z. The 
xyz coordinate system is obtained from the XY Z coordinate 
system by rotating the former around Z by $o and then ro- 
tating it around y by Qo- A typical vector R with spherical 
coordinates (R, 8, <f>) in the xyz coordinate system has polar 
coordinates (R, $) in the XY plane. 



In general, the molecules can be allowed to hop to give 
rise to Hubbard-type [25, 45] or i-J-type [46, 47] models 
with highly tunable anisotropic long-range spin-spin in- 
teractions. In this paper, however, we will assume that 
the lattice is so deep that tunneling is negligible and 
molecules are pinned in the motional ground state on 
each site. The field-free (i.e. E — 0) eigenstates of Ho are 
the simultaneous eigenstates of N 2 and N z with eigen- 
values N(N + 1) and M, respectively. Let us denote with 
\N,M) the eigenstates of H a , which the field- free eigen- 
states connect to as E is turned on. For each molecule, 
we define [28] d ± = + (<f E ± id y ), which changes the M 
of the molecule by ±1, while dP = d z couples rotational 
states with the same M. There are no selection rules on 
N for E ^ 0. 

Consider two molecules i and j separated by Ry, 
which has polar coordinates (Rij, Or,) in the XY plane 
and spherical coordinates (Rij,0ij, 4>ij) in the xyz coordi- 
nate system. The dipole-dipole interaction between these 
two molecules is [48] 

d i -d j -3(d i .R. tJ )(d j .R u ) 

H v 53 y 2 ) 

= ~W E(- 1 )' C ^9( <? «.^)2?(di > dj), (3) 

tt <-3 q=-2 

where C*{6,<f>) = yj ^p[Y kq (6 , (/)) , Y kp are spherical har- 
monics, T| 2 = dfdf, T 2 ±1 = (dldf + dfd?) A/2, and 
T 2 = (drdf + 2d° l d° 1 + d+dj) /y/6, so that T 2 changes 
the total M of the two molecules by q. T± x contributes 
only at those values of the DC electric field that have 
a ± (M — > M ±1) and ir (M-conserving) transitions of 
matching frequency. 

In each molecule, we couple states \N, M) with several 
microwave fields and choose 2S + 1 dressed states (linear 
combinations of states \N, M)) to define an effective spin 
S system. We will focus in this paper on homogeneous 



driving, which is easier to achieve as it can be done with 
microwave - as opposed to optical - fields. Extensions 
to inhomogeneous driving allow for an even richer class 
of Hamiltonians [49, 50]. Assuming dipole-dipole interac- 
tions are too weak to take the molecules out of the 25 + 1 
chosen dressed states, we can project dipole-dipole inter- 
actions onto these states to obtain a spin-S* interaction 
Hamiltonian of the form H = ^ ^2i^j H%j > where 

V(«,,.0 (/ :- II. (4) 

Here 

v(M) = (-2C 2 ,-V6Re[C 2 2 ],V6Im[C 2 2 ],Re[C 2 ],Im[C 2 ]) 

is a real five-component vector describing the five differ- 
ent angular dependences (the prefactors are chosen for 
later convenience). Each of the five components of H 
is a Hamiltonian acting on the Hilbert space of the two 
spin-S 1 particles i and j and comes with its own angu- 
lar dependence. Specific examples for the components 
of H arc given below, for example Eq. (7) for the spin- 
1/2 case. Due to Hermiticity and symmetry under the 
exchange of the two particles, each component of H has 
[(25 + l) 4 + (25 + l) 2 ]/2 independent real coefficients. 
With an appropriate choice of rotational states, DC elec- 
tric field strength, and a sufficient number of microwave 
fields, one might envision achieving full control over all 
5[(25+l) 4 +(25+l) 2 ]/2 coefficients, which, together with 
Oo, and a choice of lattice, determine the system. Re- 
quiring in addition that the total number of particles in 
any given internal state is conserved gives only (25+ l) 2 
independent coefficients in each angular dependence. 

Eq. (4) allows one to access a great variety of exotic 
spin Hamiltonians. Since we will discuss 5=1/2 exam- 
ples below in Eq. (7), here we only point out that, in Ref. 
[33], we showed how Eq. (4) can be used to realize the 
most general SU(2)-symmetric 5=1 Hamiltonian (the 
so-called bilinear-biquadratic Hamiltonian) restricted to 
the Cg angular dependence, in Rcf. [33], we also briefly 
mentioned the possibility of realizing the Kitaev honey- 
comb model. In the present paper, we discuss the details 
behind the realization of the Kitaev honeycomb model 
and provide additional insights into how to generate spin 
Hamiltonians with an arbitrary 5. 

SPIN HAMILTONIANS WITH ARBITRARY S 

In this section, we show how to obtain a variety of 
spin Hamiltonians with an arbitrary 5. For simplicity, 
in each molecule, we choose Q distinct \N,M) states 
and label them as \a), where a — 1,...,Q. We break 
this set of Q states into 25+1 disjoint sets labeled by 
p = —5, . . . , 5. We couple the states within each set with 
microwave fields to form dressed states in the rotating 
frame. We will show in detail below that ~ n microwave 
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fields are needed to couple n states and to create any 
desired linear combination out of them. We then choose 
one dressed state from each set to create the single-spin 



basis \p) = J2 a 
the 



(p) 



) with rotating frame energies 

E p . Here X) a (p) means that a is summed over the states 
belonging to the set p. The coefficients assumed to 

be nonnegative real numbers: allowing for superpositions 
with complex coefficients does not allow any additional 
tunability in this example [51]. The states \p) and \q) 
will refer to dressed states, while states \a) and |6) will 
refer to bare states. (One could also choose fewer sets 
and choose more than one dressed state from the same 
set.) For simplicity, we further assume that dipole-dipolc 
interactions are so weak and the states are chosen in such 
a way that the two-molecule state \pq) is connected via 
Hij only to itself and to \qp), while all the other pro- 
cesses are off-resonant and are negligible. Similarly, now 
in terms of the bare non-microwave-dressed states, we 
assume that the states are chosen in such a way that the 
two-molecule state \ab) is connected via only to itself 
and to 1 6a), while all the other processes are off-resonant 
and are negligible. Finally, also for simplicity, we assume 
that the states are chosen in such a way that only Cq 
contributes (an extension that includes the other four 
components of v is straightforward). The assumptions 
in the previous three sentences are generically satisfied 
if, for example, all the states \a) involved have M > 
and no accidental degeneracies occur. We then have 



1-3 cos 2 8ij 



B v\PP) (PP\ + A p A i\pi) (Pl\ 



J2^(\qp){pq\+h.c.) 



p<q 



(5) 



Where A p — J2a(p) x a^a, B p — J2a(p)<b(p) x a x bd a ,b, 

Jp,q = J2a( P )Xq) x "- x bda,b, where 

f 2^ b XM(a) = M(b) 
d a .b=l -\A,b if ma) - M(b)\ = 1 . (6) 
[ otherwise 

Here fj, a>b = (a\d M ^~ M ^ \b) and [i a = (a\d°\a). A p 
can be regarded as an effective dipole moment of state 
\p), while Bp encompasses the contribution to the inter- 
action from the transition dipole moments between the 
states a(p). On the other hand, J Piq describes the flip- 
flop transition involving states \p) and \q) and is driven 
by transition dipole moments coupling a(p) to 6(g). The 
factor of 2 and the minus sign in Eq. (6) arise from the 
fact that two identical dipoles rotating in phase in the 
x-y plane experience, on average, an interaction equal to 
negative one-half of that felt by these dipoles had they 
been pointing along z [29]. 

The implementation of SU(A r =2S'-|-l)-symmetric mod- 
els is an attractive case to consider given that such models 



have great potential for generating exotic phases includ- 
ing chiral spin liquids [52, 53]. The expression in the 
square brackets in Eq. (5) possess SU(7V=2S'+1) sym- 
metry if and only if it is proportional to ^2 p \pq) (qp\ 
up to an additive constant. Necessary and sufficient 
conditions for this are that A p A q , J p ^ q , and B p + A p 
do not depend on p and q for all p < q and that 



Bp + A p = ApA q + Jp^q/2 for all p < q. As an exam- 



ple, for N = 3, it is sufficient to satisfy A\ = A 2 = 
and B 1 = B 2 = B 3 + A\ = \.J ia = \ J lj3 = \ J 2j3 . In 
Ref. [33], we showed how to realize this and, in fact, an 
arbitrary SU(2) symmetric S = 1 interaction. Finding 
specific level configurations for higher TV is postponed 
until future work. 



THE GENERAL 5=1/2 HAMILTONIAN 

The most general S = 1/2 Hamiltonian implemcntablc 
with dipole-dipole interactions is 

RfjHij = v{6ii,<t>ij) ■ [Vmrij + W x ( ni S* + S* nj ) 
+W y { ni S] + Sfn 3 ) + W z (mS* + S? nj ) 

IT Ol QX , T QV CV I T QZ QZ 



'VV'-'i r "zz'-'i >->j 

+J xy (S?S? + SfS*) + J XZ (S*S* + sis*) 



+ 3y Z (SfSJ + SfSj)] 



(7) 



where we have 10 real 5-component vectors in the square 
brackets. By analogy with Refs. [33, 46, 47], we expect 
that these vectors might be independently controllable, 
an issue whose details we leave for future investigation. 
For this, we need ~ 50 microwave fields. In Eq. (7), rii 
is the occupation of site i. In the present manuscript, 
rii = 1 for all sites. However, this more general form of 
the Hamiltonian is necessary when hopping is allowed, 
such as in Refs. [25, 45-47]. 

In general, the interaction in Eq. (7) is also accompa- 
nied by on-site Hamiltonian terms. In particular, if the 
on-site energy difference between |^) and \l) (the two 
dressed basis states making up the spin-1/2) is much 
larger than the strength of H^, only those interaction 
terms that conserve the total S z play a role [46, 47]. In 
order to realize models such as the quantum compass 
model [54], the Kitaev honeycomb model [55], and the 
Kitaev quantum double models [56] implemented with 
two-body interactions [57], we need to break S z conser- 
vation and realize terms, such as SfSJ. To do this, we 
simply tune |f) and |4-) to be degenerate. 

For rij = 1, the Hamiltonian in Eq. (7) is the special 
case of Eq. (4): its restriction to S = 1/2. Despite the 
restriction to S = 1/2, the Hamiltonian in Eq. (7) is 
still extremely powerful, as illustrated by the following 
examples. First, the commonly used density-density in- 
teraction with the 1 — 3 cos 2 9ij angular dependence (see 
e.g. Refs. [7, 9]) comes from the first component of V, 
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while density-density interactions with more exotic an- 
gular dependences can be obtained by making use of all 
five components of V. Second, the first components of 
Jxy; Jzz, W z , and V give rise to the spin-spin, density- 
spin, and density-density interactions of the generalized 
t-J model (referred to as the t-J-V-W model) of Refs. 
[46, 47]. Third, the first three components of J xy and 
•Jzz give rise to an XXZ model featuring a direction- 
dependent spin-anisotropy [33]. In addition to provid- 
ing access to symmetry protected topological phases in 
spin ladders [33], such an XXZ model can also be used 
to obtain a variety of exotic antiferromagnets including 
a square-lattice Heisenberg model, in which the ratio be- 
tween coupling strengths on X and Y bonds is tunable 
[58]. Fourth, the 3 XX and 3 ZZ terms can be used to re- 
alize the quantum compass model [54]. In the present 
paper, we demonstrate the power of Eq. (7) on another 
important example: the Kitaev honeycomb model [55]. 



KITAEV HONEYCOMB 

In Fig. 1 and in the inset to Fig. 2(a), we show the 
Kitaev honeycomb lattice with Hamiltonian [55] 



h — J xx ' s i Sj 



-link 



y— link 



SfS] 



J. 



z — link 



We are interested in phase B, which contains the point 
Jxx — Jyy = Jzz < 0. This phase is gapless, but, in the 
presence of a magnetic field, 
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(B X S* 



ByS y 



BzS*), 



(9) 



with all three components B Xl B yi and B z being nonzero, 
acquires a gap and supports non-Abelian anyonic excita- 
tions. Non-Abelian anyons have quantum computing ap- 
plications and can be used, for example, for topologically 
protected quantum state transfer [59]. 

We use Eq. (7) for a system of n = 1 molecules per 
site on a honeycomb optical lattice - which can be im- 
plemented with three laser beams [60-62] - to implement 
the Kitaev honeycomb model. Specifically, we would like 
to obtain SfSJ interactions along <&ij = tt/6, SfSj in- 
teractions along <&ij = 7r/2, and SfS? interaction along 



It is worth noting that inter- 
I- 7r are the same. 



$ 4J = 5tt/6 [see Fig. 2(a)]. 
actions along $y and <J>y 

Let us first point out some general facts that hold inde- 
pendently of the direction and magnitude of the applied 
DC electric field and of the dressed states |f) and \\) that 
we choose. The terms in Eq. (7) involving V play no role 
since they are proportional to identity. Ignoring bound- 
ary effects, the terms in Eq. (7) involving W x , W y , and 
W 2 give rise to a uniform magnetic field, precisely as 
required by Eq. (9); we will compute this field below. 

Having taken care of the V and W Q terms, we omit 
them from the Hamiltonian. We would like the remaining 
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FIG. 2. Implementation of the Kitaev honeycomb model, (a) 
The honeycomb model: SfSj, SfSj, and SfSj interactions 
along &ij = it/6, 7i"/2, and 5n/6, respectively, (b) The rota- 
tional levels used. The diagram is schematic: the real system 
is anharmonic and levels \N,M) with the same N are non- 
degenerate (unless the levels have the same \ M\). The dressed 
states are \f) (linear combination of states indicated by ovals) 
and |4) (the rest), (c) Microwave fields that can be used to 
create the dressed state |f). Notice that level [2, 0) is used 
off-resonantly with a detuning A. (d) Microwave fields that 
can be used to create the dressed state \\). Notice that level 
1 4,0) is used off-resonantly with a detuning A. 



terms in the Hamiltonian (J xx , J yy , 3 ZZ , 3 xy , 3 XZ , 3 yz ) to 
satisfy 



Hij ({f ! f i if}) — ^3 {Sf Sj , Sf Sj , $f Sj} (10) 



with J > for the three indicated values of <Fy. Since, 
for any v($^) can be written as a linear combina- 
tion of v(7r/6), v(7t/2), and v(57r/6), Eq. (10) completely 
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determines the Hamiltonian at all ^ i 
H f$ 1 = y 1 - 2ros ( 2< ^ - 2 (" + 1 W 3 ) 

n=l d 

x fl.. ((^) (11) 

= -^3- 1 [1 - 2 cos (2* y - 47T/3)] SfSJ 
+ [l-2cos(2$ u -)]5f5j' 
+ [1 - 2 cos (2$ y - 2tt/3)] 5? 5? I . 



We will show how to realize this momentarily. First, 
observe that, for example, along $jj = 7r/3, which is 
halfway between SfS* and SfSJ, 

fly(7r/3) - -^-(2S?SJ + 2S?SJ - (12) 

Because of the l/i? 3 dependence, the strength of 
such next-nearest-neighbor interactions (along 3>jj = 
0, 7r/3, 2w/3) is reduced relative to the nearest-neighbor 
interactions by 1/3 3 / 2 w 1/5, which will be the largest 
correction introduced by long-range interactions. While 
in some cases long-range corrections are weak enough to 
ensure the survival of the desired phases [33] , it is an open 
question whether this holds for the present example. 

It now remains to find the DC electric field strength 
and direction, as well as the dressed states, that yield 
Eq. (10). We assume that the DC electric field is along 
the Z axis (i.e. O = 0) and take <J> = 0. In that case, 
the last two components of v vanish, which means that 
C±! terms do not contribute (V4 = v 5 = 0). This is fine 
because those terms are the hardest to generate exper- 
imentally as they are nonvanishing only at specific val- 
ues of the DC electric field, and it will turn out they 
are not required to generate the Kitaev Hamiltonian. 
At Go = $0 = 0, the uniform magnetic field coming 
from W Q terms is easy to calculate. Indeed, V2 and v% 
do not contribute to the magnetic field since the honey- 
comb lattice is invariant under 27r/3 rotations and since 
£LoM$i + n27r/3), « 3 ($i + n2^/3)} = {0, 0} for any 
$!. Finally, summing over all bonds connected to a given 
site, the term W • (Si + Sj) with angular dependence V\ 
gives rise to a uniform magnetic field of strength 6.58W. 

Since v 4 = v 5 = 0, we drop the last two components of 
v and H to obtain 



and 



v($) = (1, -§ cos(2$), | sin(2$)) 



H = {dy j + ±(drdf + 4dj), 



(13) 



dtdt 



d-dj 



i(dfdj-drdj)}. 



(14) 



The linear independence of v(7r/6), v(tt/2), and v(57r/6) 
allows, in principle, for the possibility of obtaining Eq. 
(10). 

We will use the 25 states shown and numbered in 
Fig. 2(b). We write the two dressed states as \a) = 
J2a(a)ya\a), where y a are complex, £ o(<t) |y a | 2 = h 
a =t,4-i and J2 a (cr) means that one sums a over the 12 
(13) states belonging to |t) (|4-)) in Fig. 2(b). We now 
keep only resonant terms and project H in Eq. (14) onto 
dressed states |f) and ||). 

We work at a DC electric field E = lOB/d. The fol- 
lowing transitions contribute to Hi. \N,M) \N',M') -> 
\N,M) \N',M'), \N,M) \N',M') \N',M') \N,M) (for 
\M — M'\ < 1, as dictated by electric-dipole selection 
rules), and \N, M) \N', —M) -> \N' , M) \N, -M). The 
transition | JV, -M- 1) \N',M) \N',-M) \N,M + 1) is 
the only one contributing to dfdj' in H2 and H^, while 
d^dj is the Hermitian conjugate of dfd^. Projecting 
these terms on |t) and |4-), we obtain an expression for 
H in terms of y a . This allows us to get an expression 
for W and for the left-hand side of Eq. (10). Since a 
magnetic field that is too large will eventually take us 
out of the Kitaev B phase, we first verify that we can 
obtain Eq. (10) with a minimal magnetic field W. One 
can then easily verify that the set of y a given in footnote 
[63] gives J = 0.0284d 2 and 6.58|W| = 7x 10~ 7 d 2 (which 
is negligibly small). 

By adjusting the coefficients y a , one can find points 
with larger values of |W|. For example, the coefficients y a 
given m footnote [64] achieve J = 0.0264d 2 and 6.58W = 
{0.00072, 0.00084, -0.03937}. 



MICROWAVE DRESSING 

So far, we have not discussed in detail what exact 
microwave couplings are required to create the dressed 
states |t) and \l) needed to implement the Kitaev hon- 
eycomb model. In this Section, we provide such a discus- 
sion. The methods we present here are generally appli- 
cable to the creation of generic spin models from dipolar 
interactions. 

Although generating large numbers of precisely tuned 
microwaves is experimentally feasible [65], the effort in 
obtaining and verifying the spectrum of microwaves in- 
creases with the number of frequencies. It is thus ex- 
perimentally desirable to use as few microwave fields as 
possible. In order to create a dressed state out of n bare 
states, one needs at least n — 1 microwaves (assuming the 
same microwave cannot be used to couple more than one 
transition). Therefore, let us first address the question 
of whether n — 1 microwaves suffice to prepare a dressed 
state at any prescribed rotating frame energy featuring 
any prescribed superposition of the bare states. 

Suppose we would like to prepare a dressed state \D) = 
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S"=i V a \ a ) ou ^ °f the n bare states \a) (a = 1, . . . ,n) 
with prescribed coefficients y a and a prescribed rotating- 
frame energy E relative to the bare state |m). Sup- 
pose state | a) is coupled to state |a + l) using a mi- 
crowave with Rabi frequency fl a for a = 1, ... ,n — 1. 
The rotating-frame Hamiltonian is then 



/Ai fti 
n$ A 2 2 




n* 2 a 3 











A„-i 






a„; 



2/2 

2/3 



2/n-l 



2/2 

2/3 



.(15) 



Here A m = (recall that \m) is the reference sate), while 
A a+ i — A a is the detuning of the a'th microwave field 
from the \a) — > \a + 1) transition. Writing 

O a = fZ a exp{iarg(y Q /y a+ i)} , (16) 

where f2 a are real, we obtain n real-number equations 
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2/n-l 












Vn 



E - A 2 , 

£-A„, (17) 

S-A„_i, 
£7- A„. 



One can satisfy the equations using n— 1 Rabi frequencies 
fi a and any one nonzero detuning A c . Allowing other 
detunings to be nonzero, one obtains freedom in how to 
choose the magnitudes of il a . This can be useful, for 
example, in cases when the magnitudes of fl a need to be 
larger than a certain value (e.g. larger than the strength 
of hyperfine interactions, as in Refs. [33, 47] - see our 
Outlook). 

The approach we have just described works nicely in 
the cases where the only off-diagonal bare-states pro- 
cesses that contribute to the final spin Hamiltonian have 
the form \a) \b) — > \b) \a). This holds for all the config- 
urations studied in Refs. [46, 47]. It also holds [66] for 
the level configuration used to implement the bilinear- 
biquadratic model in Ref. [33]. 

On the other hand, consider a model where , 14-)} = 
{|0,0),j/_i |1,-1) +yi |1,1)} (as in Ref [33]) and where 
we want the process |1,— 1)|0, 0) — > |0, 0) j 1, 1) to con- 
tribute. We can couple states [1,— 1) and [1, 1) with two 
microwave fields via state |2,0). It is crucial that these 



two microwave fields have the same frequency; other- 
wise, in the rotating frame, the process |1,— 1) |0, 0) — > 
1 0, 0) |1, 1) would pick up an oscillatory time dependence 
and will average out to zero. This means that, in the 
cases where off-diagonal bare-state processes other than 
| a) | b) — > \b) | a) contribute, there are restrictions on what 
detunings can be chosen as nonzero. 

In the case of the configuration for the Kitaev honey- 
comb model shown in Fig. 2(b), it is easy to check that 
the microwave couplings shown in Fig. 2(c) and Fig. 2(d) 
for |t) and [4-), respectively, work. In particular, in Fig. 
2(c) all microwaves are assumed to be resonant except 
for the microwaves coupled to state |2,0), which are all 
detuned by the same amount A. This detuning provides 
the extra tuning parameter one has to add to the Rabi 
frequencies in order to fully control the coefficients i/j 
defining state |f)- The fact that [2, 0) has M — ensures 
that it is nondegenerate. Hence, any contributing bare- 
state process containing |2, 0) in the initial two-molecule 
state will also contain it in the final two-molecule state, 
meaning that a detuning of this state keeps the process 
resonant. Similar arguments apply to Fig. 2(d), where 
all microwaves are assumed to be resonant except for the 
microwave coupling to state |4, 0). To ensure that |t) and 
\l) are degenerate, we use [2, 2) and |2, —2) as the two 
reference states and set the two drcssed-state energies to 
be equal. It is also worth pointing out that an additional 
N = 5 state is introduced in both Fig. 2(c) and Fig. 2(d) 
to act as an intermediate state for coupling two states 
with AM = 2. We also note that one has to be careful 
using 7r transitions to couple states with M/0 since the 
same microwave will act on both M and — M. 

A solution alternative to the introduction of a detun- 
ing is an introduction of additional resonant microwave 
couplings, resulting in configurations such as three states 
coupled with three microwave fields. At the same time, it 
is worth pointing out that, in some cases, one may wish 
to introduce additional detunings on purpose with the 
goal of eliminating certain interaction processes (writ- 
ten in terms of bare non-microwave-dressed states as 
|a, b) — > \a',b')), which can be regarded as an additional 
control knob. 



OUTLOOK 

The most straightforward implementation of our pro- 
posal would involve polar molecules with no nuclear spin. 
This is possible, for example, in the case of a SrYb 
molecule, which seems to be within reach given the avail- 
ability of quantum degenerate Sr and Yb gases [39]. In 
the case where hyperfine structure is present, such as in 
alkali dimers, we assume that the Rabi frequencies f2 a arc 
larger than the strength of hyperfine interactions H h{ . In 
situations where this poses a significant restriction, the 
opposite limit, in which flj <C -Hhf, can also be consid- 
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ered. 

While we have focused on pinned molecules, allowing 
the molecules to hop in the lattice gives rise to t-J-type 
[46, 47] or Hubbard-type [25] models with highly tunable 
long-range spin-spin interactions featuring a direction- 
dependent spin anisotropy. As a particularly simple ex- 
ample, one can obtain a single-component system with 
a density-density interaction whose angular dependence 
is given by a linear combination of all five Cg [see the 
V term in Eq. (7)]. This may allow one to study, for 
example, exotic quantum Hall physics [67]. 

The approaches we presented can be used to tune 
dipole-dipole interactions to zero [12, 29]. In that case, 
one has to consider corrections due to dipole-dipole inter- 
actions to second order (which give van-der-Waals-type 
1 /R 6 dependence) and interactions between higher multi- 
pole moments [68] . These interactions can also be viewed 
as extra control knobs. 

We have focused in the present manuscript on polar 
molecules. However, the ideas that we presented should 
be extendable in a straightforward manner to other sys- 
tems interacting via dipole-dipole interactions such as 
magnetic atoms [69, 70], Rydberg atoms [71, 72], and 
magnetic solid-state defects [73, 74]. 
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